Periodic Pattern Formation in Evaporating Drops 
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Solute deposits from evaporating drops with pinned contact line are usually concentrated near 
the contact line. The stain, or pattern, left on the substrate then consists of a single ring, commonly 
known as a coffee ring. Here we report on a variation of this phenomenon when periodic patterns 
emerge. We attribute these to phase transitions in certain solutes as solute concentration increases. 
Examples may include dissolved to crystalline transition in salt or order- disorder transition in liq- 
uid crystals. Activated nature of the phase transitions, along with the newly imposed boundaries 
between phases, may then invert solute density profile and lead to periodic deposits. Hereby we 
develop a general theoretical model and report on experimental observations on salt in water. 
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Evaporation rate of sessile drops is generally limited by 
the rate of vapor diffusion in the surrounding atmosphere 
P, . Faster evaporation near the contact line produces 
outward radial flow of liquid and solute. Residual solute 
then accumulates near the contact line and a ring of so- 
lute deposit is left on a substrate. This is known as the 
coffee-ring effect [H, H, 0] ■ Here we argue that situation 
changes once solute, such as salt or liquid crystal, allows 
for phase transitions. Once solute concentration exceeds 
the saturation density, solute suspension becomes unsta- 
ble and a new phase is nucleated. Diffusion of solute near 
the phase boundary becomes relevant, along with the liq- 
uid induced solute flow, and concentric deposit rings may 
emerge. 

Formation of solute rings may clearly be seen in salt 
deposits left by evaporating salt-water mixtures (Figure 
[I}. For this experiment ternary solutions of sodium chlo- 
ride (NaCl) and sodium polystyrenesulfonate (NaPSS) 
in deionized water were used [f|. Fixed- volume droplets 
were deposited onto glass substrate and allowed to evap- 
orate. Salt deposits were then visualized with an optical 
microscope. We attribute the role of NaPSS in these 
experiments to pinning the droplet contact line, as dis- 
cussed later in the text. 

Phase transitions render homogeneous solute suspen- 
sion unstable above the critical concentration. Crystal 
nucleation or spinodal decomposition result in formation 
of new phases, such as salt crystals of Figure [TJ In the 
evaporating drops, the first of these phase transitions 
takes place near the rim of the drop, in the region of 
the highest solute concentration. This is similar to the 
coffee-ring effect [H, H, H| . Similarity, however, ends here 
as further flow of solute is affected by the boundary be- 
tween phases. On the one hand, higher evaporation rate 
and evaporation induced flow of solute support higher 
solute concentration near the contact line. On the other 
hand, solute concentration at the phase boundary may 
not exceed the critical concentration. The solute con- 
centration profile may then invert, with concentration 
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FIG. 1: Concentric rings observed in salt deposits of drying 
droplets. Drop diameter is approximately 2 mm. 



decreasing near the contact line. This is schematically 
illustrated in Fig. [2^. Location of the peak in the solute 
concentration is the likely place for the subsequent ring 
nucleation. 

Proposed mechanism implies existence of the oversatu- 
ration region. Higher concentration of solute, often called 
supercritical concentration @, 0|, is needed for reliable 
phase nucleation. At concentrations just above critical, 
nucleation is possible, yet very slow. At the supercritical 
concentration, on the other hand, nucleation times be- 
come relevant relative to the evaporation times. Solute 
concentration may even reach the spinodal point so that 
phase transitions occur spontaneously. We will use the 
term nucleation irrespective of the actual origin of the 
supercritical concentration. 

We now qualitatively estimate ring spacing Ar. Dif- 
fusive flow of solute near the contact line is DAtp/Ar, 
where A<p is difference between critical and supercriti- 
cal concentrations. Evaporation induced flow is uo<p ~ 
4>R/Tq, where uq = R/Tq is characteristic velocity of 
the liquid, and R and To are radius of the drop and 
its evaporation time. Rings emerge when two flows be- 
come comparable, so that the typical ring spacing is 
Ar ~ DTo A<fi/R(j). This is of the order of tens of microns 
for salt in water. Concentric rings of this separation are 
readily resolved (Fig. [T]). 

Pattern formation described here complements other 
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FIG. 2: Schematic illustration of the concentric ring formation. Different curves represent solute density profiles at different 
times, with higher curves corresponding to progressively later time steps. Critical solute concentration for phase transitions 
is set at 1.0, and initial concentration is set at 0.75. Image a) illustrates increase in the solute concentration next to the first 
crystalline ring, assumed to be at radial coordinate f = 0.95. Solute concentration is fixed at critical value near the crystalline 
ring. Elsewhere in the drop concentration progressively increases with time producing oversaturated region; b) Once the peak 
concentration reaches suercritical value (<f> c — 1.20 in this example), new crystalline ring is formed; at this moment phase 
boundary is moved to the location of the newly formed ring; c) Formation of progressive crystalline rings in the drop. 



non-equilibrium phenomena observed in evaporating 
drops. Convective flows were speculated to modulate 
polymer deposits in evaporating drops [8[ . Surface forces 
and surface rupture were reported to order submicron 
sized particles [l, [Io|, [H| and salt deposits [l2[ in drying 
drops. Viscous drag was found to align DNA molecules 
near the contact line [l3[. Yet none of these mecha- 
nisms deals with phase transitions in solutes. Proposed 
mechanism is also distinctly different from the Liesegang 
ring formation in gels [f| 0] , which is also based on salt 
crystal nucleations. In the evaporating drops nucleation 
events do not require diffusion of two solute species, as 
in Liesegang rings. Instead, concentration gradients and 
nucleation events arise from the competition between so- 
lute diffusion and evaporation-driven liquid flow. 

We now formulate our model. Comparison with exper- 
imental observations will follow. Sessile drop evaporation 
rate is primarily determined by the diffusion rate of liq- 
uid molecules in the atmosphere surrounding the drop, 
rather than by the rate at which these molecules escape 
the surface of the drop. This leads to the fundamental 
prediction that evaporation rate is highest next to the 
contact line, and decreases as one approaches the cen- 
ter of the drop l]. In the limit of small contact angles 
and low solute concentrations the evaporation rate J(r) 
is given by [3| 
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where D a i r is diffusion constant for vapor particles in 
the atmosphere surrounding the drop, p v and are 
saturated and ambient vapor densities, and R and r are 
drop radius and distance from the center respectively. 
Throughout this paper we assume that concentration of 
solute is small enough not to affect liquid flow or evapo- 
ration rate. 

The total evaporation time for a drop with initial con- 
tact angle 0q and liquid density p is 



Tn = 
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This time is much longer than typical relaxation times 
inside the liquid. The drop maintains its spherical cap 
shape throughout the evaporation, with height profile 
h(r) — 8(R 2 — r 2 )/2R. Here 9 is contact angle, assumed 
small. 

Radial velocity u of the liquid in the drop may be read- 
ily evaluated from the constitutive equation 



-V(rhu) + - + d t h = 0, 
r p 

which immediately integrates to 
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So far we have only considered flow of liquid. Solute 
particle, that generally move with the liquid, may also 
diffuse in it with concentration profile (f> satisfying 
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where D is diffusion coefficient for solute particles. For 
the time being we assume that flow is one dimensional 
and ignore azimuthal coordinates. It is convenient to 
introduce dimcnsionlcss time t = t/To and coordinate 
f = r/R. The main equation for the solute profile is 
therefore 



dt [(1 - t)cf\ = 
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Typical value for the parameter T = D ( P ' - P — ) p o 
room temperature for salt diffusion in water is T ~ 0.7#o, 
assuming p w 10 3 kg/m 3 , p v — p^, ps 1.7 • 10~ 2 kg/m 3 , 
D air w 0.24 • 1CT 4 m 2 /s, D rs 1.5 • 1CT 9 m 2 /s. 

At the onset of evaporation, solute flow follows the 
well-known coffee- ring phenomenon 0]. Diverging evap- 
oration rate at the contact line creates excess in solute 
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FIG. 3: Numerical calculations of ring spacing Af. a) Change 
in ring spacing with F; curves, top to bottom, correspond to 
4> c = 1.2 and 4>o = 0.95; (f> c = 1.1 and 0o = 0.95; C = 1.1 and 
4>o = 0.75. b) Change in ring spacing with <j> c . Curves, top 
to bottom, correspond to 0o = 0.95 and V = 0.03; <f)o = 0.95 
and T = 0.02; <f) = 0.75 and T = 0.02. 



concentration, and the first ring nucleates near the con- 
tact line. We do not consider this effect here. Further 
flow of the solute proceeds with a new boundary condi- 
tion (j)\f.—i_s = 1, which requires that solute concentra- 
tion near the newly formed dense phase must equal the 
critical concentration. Here 8 is location of the nucleated 
ring, and we have rescaled all solute concentrations in 
units of the critical concentration. 

The most interesting aspect of the deposit is the spac- 
ing between rings. We start with numerical analysis, and 
calculations are carried out as follows. The initial density 
concentration is set at certain value between <fia < 1, 
and adsorbing boundary condition is placed at f = 1 — 5. 
The system is then allowed to evolve according to Eq. [6] 
until concentration reaches supercritical value <p c at some 
radial coordinate f (Fig. ^ . Next crystalline ring is pre- 
sumed to form at this position and adsorbing boundary 
is moved to the location of the newly formed ring. Cal- 
culations are then allowed to proceed. This process is 
schematically illustrated in Fig. O 

In the numerical analysis, we varied parameters in 
the range T = 0.002 - 0.05, <f> = 0.05 - 0.95, and 
(f> c = 1.05 — 1.40. The grid size was set small enough 
to get stable solution, typically around 10 -4 for coordi- 
nate step, and 10 -8 for time step. The ring spacing Af, 
which was found surprisingly uniform across large area 
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FIG. 4: Comparison of numerically calculated ring spacing 
(vertical axis) versus predictions of interpolating expression 
((9J (horizontal axis). 



of the drop surface, was measured as a function of these 
parameters, as shown in Fig. [3] Particular attention was 
given to the range of spacings Af ~ 0.01 — 0.05, which 
corresponds to experimentally achievable spacing for salt 
deposits discussed later in the text. Outside this range 
experimental system becomes unstable towards periodic 
pattern formation. 

The first conclusion is that predicted ring spacing is 
roughly linear in T. This scaling can actually be ex- 
tracted from the asymptotic expansion of Eq. O We first 
introduce new variable z = \/\ — r 2 , which is closely re- 
lated to the distance from contact line, and expand Q 
in small z 
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This equation can be converted into parameter free form 
using substitution z. = Tz and t — T 3 t. Leaving terms 
small in time and T, 
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Spacing between rings at fixed location f in the drop 
should therefore scale as Af ps Azy/1 — f 2 /f oc T, in line 
with the numerical calculations (Fig. [3^l). 

Concentration dependence of the ring spacing is richer. 
Experiments and numerical calculations show that ring 
nucleation events proceed in rapid succession. At that 
time average solute concentration is actually close to crit- 
ical concentration (1.0 in our rescaled variables). That 
would imply that 1 — f ~ 1 during early ring nucle- 
ation, and spacing should therefore scale as Af oc TtpQ. 
Fitting to numerical data shows that this scaling is in- 
deed close to the actual trend, with Af ~ A((j>o)T, where 
A((f> ) w cj> - 0.036 - 13.97(0 C - l)r. 

Finally, we can add dependence on the supercritical 
concentration <p c . Assuming that ring spacing varies as 



Af = A{(j> )B(cp c 



nr 



(9) 



we find B{(j) c - 1) w 0.874[1 + 14.42(0 C - 1) - 19.1O(0 C - 
1) ]. Comparison between this interpolating expression 
and full numerical calculation is given in Fig. 0J 
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FIG. 5: Ring spacing in evaporated salt deposits. Open sym- 
bols correspond to NaPSS concentration 1 g/L, and filled sym- 
bols correspond to 0.1 g/L. The shape of each symbol is as- 
sociated with the salt concentration, as shown in the legend. 
Uncertainty in measured values is estimated at 10%. Solid 
line is a linear fit to the data points. 

Aforementioned predictions may now be compared 
with the experiments on salt deposit growth in evapo- 
rating salt-water mixtures. For this purpose we have 
used ternary solutions of sodium chloride (NaCl) and 
sodium polystyrenesulfonate (NaPSS) in deionized wa- 
ter. Fixed- volume droplets of this solution were de- 
posited onto glass substrate and allowed to evaporate in 
a humidity-controlled glovebox. Glass substrates were 
a priori rinsed with RO water for a few minutes and 
dried in nitrogen atmosphere. Observed patterns (Fig. [1]) 
where then analyzed for the ring spacing. We attribute 
the role of NaPSS in these experiments to pinning the 
droplet contact line. Neither degree of polymerization 
of NaPSS (varied in the range between 100,000 Da and 
1,000,000 Da) nor its concentration (0.1 g/1 to 1 g/1), 
seemed to affect the ring spacing. We refer the reader to 
Ref. H for detailed discussion of the experimental proce- 
dures. 

According to the predictions above, the spacing be- 
tween salt rings is expected to follow Ar oc T4>qR oc 
Oq^qR oc vo(f>o/R 2 , where vq, 4>q, and R are initial vol- 
ume of the droplet, salt concentration, and radius of the 
droplet respectively. To test this prediction, we used 
drops with initial volume 0.5 fiL and varied initial salt 
concentration between 120 and 300 mM. Additionally, 
glass substrates of varying cleanness were used to allow 
for larger spread in drop radii. Observed drops had radii 
in the range 1.0 - 1.3 mm. The results are summarized 
in Fig. and appear in good agreement with the above 
predictions. The uncertainty in the measured ring spac- 
ing and drop radius is estimated at 10%. 

Observed trends in salt ring spacing may not be read- 



ily explained by other scenarios of ring formation. For 
example, nucleation of salt rings could cause surface rup- 
ture and contact line depinning (T^. Dissolved polymer 
in the drop may then create repetitive pinning effect. In 
such a scenario, however, ring spacing would likely de- 
pend on polymer concentration rather than on salt con- 
centration, contrary to the observations. Also, surface 
rupture is likely to favor unstable, fractal-like patterns 
pdl . [l5| . Another scenario could come from instabilities 
such as convective flow or Marangoni flow ||. Yet these 
flows are likely to be suppressed near the contact line, 
where polymer concentration is very large. These flows 
are also likely to be sensitive to polymer concentration, 
contrary to the observations @ . 

Finally, a few words on the stability of the concentric 
rings. In the discussion so far we have assumed radial 
symmetry of the problem. Yet nucleation of the crys- 
talline rings occurs in the oversaturated regions of the 
drop. Growth of new phase should then be characterized 
by fingering instabilities of dendritic growth [lj|, [I?], 03 • 
In evaporating drops these instabilities are indeed possi- 
ble, yet their radial extent is suppressed by the width of 
the oversaturated region (see Fig. ^jp). Rapid decrease 
in solute concentration away from the supercritical point 
prevents instabilities from radial growth. 

Situation changes at late stages of evaporation, when 
the whole volume of the drop becomes oversaturated 
(Fig. [2k) . Fingering instabilities are no longer con- 
strained and may extend through the whole drop. This 
is indeed observed in salt deposits: in the inner region 
of the drop periodic ring formation is clearly replaced by 
radially oriented salt crystals or other disordered struc- 
tures (Fig. [l}. The transition between periodic rings and 
irregular structures in the inner region is common to all 
drops considered. 

In conclusion, we have shown that phase transitions in 
solute may qualitatively change deposit growth in evap- 
orating drops. The phase equilibrium between solute 
phases may invert the solute concentration profile, and 
lead to the growth of nearly periodic patterns. Formation 
of a single coffce-ring deposit is replaced by a sequence 
of concentric rings spanning large part of the drop area. 
The spacing between rings may be controlled via changes 
in ambient atmosphere, temperature, diffusion constant, 
and other parameters. This may prove crucial in under- 
standing many contact line phenomena in evaporating 
liquids, and lead to better control over pattern growth. 
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